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We derive an equation determining the small- a; evolution of the F2 structure function of 
a large nucleus which includes all multiple pomeron exchanges in the leading logarithmic 
approximation using Mueller's dipole model jl|-0|. We show that in the double leading 
logarithmic limit this evolution equation reduces to the GLR equation . 
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I. INTRODUCTION 



The problem of understanding the large gluon density regime in high energy scattering has always been 
one of the challenges of perturbative QCD (PQCD). Unitarity of the total cross-section and saturation of the 
gluon distribution are among the most important issues related to the problem. The BFKL equation 
is the only well-established tool of PQCD which allows us to explore this high density region by resumming 
the leading longitudinal logarithmic contribution to the scattering process. In BFKL evolution the small- 
x partons are produced overlapping each other in the transverse coordinate space ||, therefore creating 
1-^ ' high density regions in the hadron's wave function (hot spots). The next-to-leading order correction to 

BFKL equation has been calculated recently [ ^0|JTl|l . Although the final conclusion one should draw from 
the calculation of 10 is still to be understood, there are some serious problems associated with the 



interpretation of the result |T^-|l^|. However, we are not going to address these issues in this paper for the 
following reason. As was shown in jl2|,[l3| the effects of the second order BFKL kernel become important 
in hadron-hadron scattering at the rapidities of the order of Y/vlo ~ 1/a 5 / 3 , with a the strong coupling 
constant. At the same time the unitarity constraints, associated with the multiple (leading order) hard 
pomeron exchanges are expected to be reached at Yg ~ (l/a) ln(l/a) |lq ], which is parametrically smaller 
than Yjsslo for small coupling constant. Therefore multiple pomeron exchanges become important at lower 
center of mass energies than the effects of subleading corrections, possibly leading to unitarization of the 
total hadron-hadron cross-section. Hence the problem of resummation of the multiple pomeron exchanges 
seems to be more important for describing the recent experimental results, such as ZEUS 1995 data Eg ], 
which probably shows evidence of saturation of the F2 structure function at low Q 2 . 

In this paper we are going to consider deep inelastic scattering (DIS) of a virtual photon on a large nucleus 
and will resum all multiple pomeron exchanges contributing to the F2 structure function of the nucleus in 
the leading longitudinal logarithmic approximation in the large iV c limit. The first step in that direction 
in PQCD was the derivation by Gribov, Levin and Ryskin (GLR) of an equation describing the fusion of 
two pomeron ladders into one in the double logarithmic approximation || . The resulting equation with the 
coefficient in front of the quadratic term calculated by Mueller and Qiu |6| for a low density picture of a 
spherical proton of radius R reads 



d\n(l/x)d\n(Q 2 /A 2 QCD ) TT x ^ ' 3C F R 2 Q 2f 
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This equation sums up all multiple hard pomeron exchanges in the gluon distribution function in the double 
logarithmic limit. 

Since then there have been several attempts to generalize the GLR equation. Recently an equation has 
been proposed by Ayala, Gay Ducati, and Levin in [jl7|,[l8|, which tries to incorporate the Glauber-type 
multiple rescatterings of a probe on the nucleons in a nucleus (see Fig. Q). Using the results of Mueller in 
p9[ for a pair of gluons multiply rescattering inside a nucleus, the authors of Jl7],[l8| obtained the following 
equation for the gluon distribution of the nucleus in the double logarithmic approximation 



d 2 xG A {x,Q 2 



N C C F S. 



d ln(Q 2 / h? QCD )d ln(l/ x) 



exp 



2ait 2 1 



xG A (x,Q 2 ) 



(2) 



If one expands the right hand side of Eq. (|2j) to the second order in xGa one recovers the GLR equation 
[Eq. (^) is written here for a cylindrical nucleus. Therefore the coefficients in the obtained GLR equation 
will not match those of Eq. d3)]. 
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FIG. 1. Forward amplitude of DIS on a nucleus in the quasi-classical (Glauber) approximation. 



An extensive work on resumming the multiple pomeron exchanges in the gluon distribution function in 
the leading ln(l/a;) approximation (i.e. without taking the double logarithmic limit) has been pursued by 
Jalilian-Marian, Kovner, Leonidov, McLerran, Venugopalan, and Weigert [20|-|2^1. Starting with a model of 
a large nucleus |2(||26|], which provides some effective action [^0 21 1, they develop a renormalization group 
procedure which integrates out harder longitudinal gluonic degrees of freedom in the nucleus and allows one 
to resum leading ln(l/x) contribution to the gluon distribution function. The resulting equation, written 
in a functional form in p3fl is supposed to resum these leading logarithms including also all powers of color 
charge density of the nucleus, which in a more traditional language corresponds to resummation of multiple 
pomeron exchanges. However, even though that equation at the lowest (one pomeron) level reduces to the 
expected BFKL equation |pffi| , in general it is a very complicated functional differential equation which can 
not be solved even numerically. Recently the double logarithmic limit of that equation was obtained pj| ], 
providing us with another equation for xG: 



d 2 xG(x,Q 2 ) 



N C (N C - 1) 



dHQ 2 /A 2 QCD )dHl/x) 



S±Q 2 



1 n(N c -l)Q 2 S ± 
1 + 2^G ° XP 



tt(N c -1)Q 2 S^ 
2axG 



xEi 



ir(N c - 1)Q 2 S ± 
2axG 



(3) 



Eq. (3) is written here for a cylindrical nucleon with transverse area Sj_. In the limit of small gluon density 
Eq. (3) reduces to GLR equation ]25| . 

Our approach in this paper is pursuing the same goals as authors of ]20|-|25[] . We will also write an 
equation which resums all multiple pomeron exchanges on a nucleus in the leading logarithmic approximation. 
However our strategy is a bit different from pOl-pl . We will consider the scattering of a virtual photon on 
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a nucleus at rest, therefore putting all the QCD evolution in the wave function of the virtual photon. This 
is different from what was done by the authors of p0|-p5|, since they were developing the QCD evolution 
inside the nucleus. The virtual photon's wave function including the leading logarithmic evolution was 
constructed in the large iV c limit by Mueller in . This so-called dipole wave function in fact contains all 
multiple pomeron exchanges, which in the large N c language can be pictured as multiple color "cylinders" . 
A numerical analysis of the unitarization of the total onium-onium cross-section through multiple pomeron 
exchanges was carried out in the framework of the dipole model by Salam in . Considering the scattering 
of a virtual photon (quarkonium state) on a nucleus simplifies the problem, allowing us to treat it analytically. 
In Sect. II we will use the dipole wave function to write down an equation which governs the evolution of the 
i*2 structure function of the nucleus (formula Jl5|)). Our equation is directly related to a physical observable 
(F 2 ) and, therefore, is free from all the problems and ambiguities associated with dealing with the gluon 
distribution function xG. It is a non-linear integral equation, not a functional differential equation like in 
p2| . Therefore one should be able to solve Eq. ( |l5| ) at least numerically. 

We will dedicate Sect. Ill to exploring the equation resulting from taking double logarithmic (large Q 2 ) 
limit of the equation derived in Sect. II. We see that our equation reduces to GLR equation, failing to 
reproduce Eqs. (||), (||). Finally, in Sect. IV we will conclude by discussing the limitations of the large N c 
approximation, as well as some advantages of our approach. 



II. EVOLUTION EQUATION FROM THE DIPOLE MODEL 

We start by considering a deep inelastic scattering process on a nucleus. As shown in Fig.|l| , the incoming 
virtual photon with a large q+ component of the momentum splits into a quark-antiquark pair which 
then interacts with the nucleus at rest. We model the interaction by no more than two gluon exchanges 
between each of the nucleons and the quark-antiquark pair. This is done in the spirit of the quasi-classical 
approximation used previously in |is|]28| , |29fl . The interactions are taken in the eikonal approximation. Then, 
as could be shown in general, i.e., including the leading logarithmic QCD evolution, the total cross-section, 
and, therefore, the F 2 structure function of the nucleus can be rewritten as a product of the square of 
the virtual photon's wave function and the propagator of the quark-antiquark pair through the nucleus 
@0 .The expression reads Q 

F 2 (x,Q 2 ) = —^ / — ^-[$ r (xoi,*) + 4> L (x 01 , z)] d 2 b AT(x 01 , b , Y), (4) 

where the incoming photon with virtuality Q splits into a quark-antiquark pair with the transverse coordi- 
nates of the quark and antiquark being xo and xi correspondingly, such that xio = xi — xo. The coordinate 
of the center of the pair is given by bo = ^(xi + xo). Y is the rapidity variable Y = \ns/Q 2 = lnl/x. 
The square of the light cone wave function of qq fluctuations of a virtual photon is denoted by $t(xoi, z) 
and <!>l(x:oi,z) for transverse and longitudinal photons correspondingly, with z being the fraction of the 
photon's longitudinal momentum carried by the quark. At the lowest order in electromagnetic coupling 
(01 em) &t( x oi,z) and $i(x i,z) are given by [ pO| , ^l| and references therein] 

$ T (x 01 ,z) = | a 2 K*(x 01 a) [z 2 + (1 - z) 2 ]} , (5a) 



<Mx 01 ,z) = 2 J^EI± 4 QV (1 - zf K 2 (x 01 a), 



(5b) 



with a 2 — Q 2 z(l — z). We consider massless quarks having only one flavor. 

The quantity Njxpi, bp, Y) has the meaning of the forward scattering amplitude of the quark-antiquark 
pair on a nucleus |31|. At the lowest (classical) order not including the QCD evolution in rapidity it is given 
by 



A(x O i,b o ,0) = -7(x i,b ) 



1 — exp 



2X 



(6a) 
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with v as defined in |28(] and A being the mean free path of a gluon in a nuclear medium, as defined in |2S 
In the logarithmic approximation for large Q 2 (small xoi) Eq. (|6a|) can be rewritten as 



iY(x i,b o ,0) = -7(x i,b ) ~ < 1 - exp 



2N C S±_ 



X^ x AxG{x, 1/x; 



(6b) 



7(xoi,bo) is the propagator of the qq pair through the nucleus. The propagator could be easily calculated, 
similarly to |^8|,^0|, giving the Glauber multiple rescattering formula (||). Here and throughout the paper we 
assume for simplicity that the nucleus is a cylinder, which appears as a circle of radius R in the transverse 
direction and has a constant length 2R along the longitudinal z direction. Therefore its transverse cross- 
sectional area is S± — irR 2 . In formula (JsJ) A is the atomic number of the nucleus, a is the strong coupling 
constant and xG{x, 1/x^) is the gluon distribution in a nucleon in the nucleus, taken at the lowest order in 
a, similarly to Pq] , 

Eq. (^) resums all Glauber type multiple rescatterings of a qq pair on a nucleus. As was mentioned before, 
since each interaction of the pair with a nucleon in the nucleus is restricted to the two gluon exchange, the 
formula (^) effectively sums up all the powers of the parameter a^A 1 / 3 . Or, looking at the power of the 
exponent in (^) we conclude that since xqi ~ l/Q, it resums all the powers of a q-i ■ This is the definition 
of quasi-classical limit, a more detailed discussion of which could be found in p6| . 

Since the nucleus is at rest in order to include the QCD evolution of F2 structure function, we have to 
develop the soft gluon wave function of the incoming virtual photon. In the leading longitudinal logarithmic 
approximation (In 1/x) the evolution of the wave function is realized through successive emissions of small-x 
gluons. The qq pair develops a cascade of gluons, which then scatter on the nucleus. In order to describe 
the soft gluon cascade we will take the limit of a large number of colors, N c — > 00. Then, this leading 
logarithmic soft gluon wavefunction will become equivalent to the dipole wave function, introduced by 
Mueller in [|]-[|. The physical picture becomes straightforward. The qq pair develops a system of dipoles 
(dipole wave function), and each of the dipoles independently scatters on the nucleus, as shown in Fig. ^. 
Since the nucleus is large, we may approximate the interaction of a dipole (quark-antiquark pair) with the 
nucleus by j(x, b) given by Eq. (Q), with x and b being the dipole's transverse separation and impact 
parameter. That means that each of the dipoles interacts with several nucleons (Glauber rescattering) in 
the nucleus independent of other dipoles. The interaction of each of the dipoles with the nucleus is the same 
as was shown in Fig. [l] for the initial qq pair. 




FIG. 2. DIS on a nucleus including the QCD evolution in the quark-antiquark pair in the form of dipole wave 
function. Each double line represents a gluon in the large-iV c limit. 

To construct the dipole wave function we will heavily rely on the techniques developed in [[l] |4| . Following 
JjL|| we define the generating functional for dipoles Z(bo, X01, Y, u) (see formulae (16) and (17) in 0). The 
generating functional then obeys the equation (see Eq. (12) in f|) 
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Z(bo,x i, Y,u) = u(b ,x i)exp 



(iy exp 



d 2 x 2 



' 01 



2 2 
X 02 X 12 



Z(b + ixi2,x 2,y,'u) Z(b - ix20,Xi2,J/,«), 



(7) 



where X20 = xo — x 2 , X21 = Xi — X2 and the integration over X2 is performed over the region where X02 > p 
and X12 > p. This p serves as an ultraviolet cutoff in the equation and disappears in the physical quantities, 
bo = ^(xo + xi) is the position of the center of the initial dipole in the transverse plane M. Cf — N c /2 in 
the large N c limit. The generating functional is defined such that Z(bo,Xoi,y, u = 1) = 1 [see jl]]]. 



Analogous to we now define the dipole number density by 



2irx' 



-ni{x m ,Y, |b - b |,a:) = 



(5u(b, x 



■Z(bo,x i,y,u)| u= i 



(8) 



n\{xQ\, Y, |b — bo|,a;) convoluted with the virtual photon's wave function gives the number of dipoles of 
transverse size x at the impact parameter |b — bo| with the smallest light cone momentum in the pair 
greater or equal to e~ Y q + . Similarly to the dipole number density we can introduce dipole pair density ||.|| 
for a pair of dipoles of sizes x\ and X2 at the impact parameters |bi — t>o | and|b2 — t>o | by 



1 



1 



2itx\ 2ttx. 



rn 2 (x i,Y, |bi - b |,a;i, |b 2 - bo|,ai2) = 



2! <$u(bi,xi) <5u(b 2 ,x 2 ) 



Z(b ,x i, Y", u)\ u=1 . 



(9) 



Our notation is different from the conventional approach of J@|| by the factor of a factorial, for reasons 
which will become obvious later. Generalizing the definition (M) to k dipoles of sizes x\, . . . ,x n situated at 
the impact parameters | t>i — bo|, . . . , |b& — bg| we easily obtain: 



II ~^~2 nk ( Xol > Y > l bl ~ b o\,Xi, |b fe - b \,x k ) = ^_ Y[ 



Z(b ,x 01 , Y~, u)\ u=1 



(10) 



One can now see that in order to include all the multiple pomeron exchanges one has to sum up the 
contributions of different numbers of dipoles interacting with the nucleus. Namely we should take the dipole 
number density rti(a;oi, Y, b, x) and convolute it with the propagator of this one dipole in the nucleus 7(x, b). 
Then we should take the dipole pair density n2(xoi, Y, bi, x\, \>2, £2) and convolute it with two propagators 
7(xi, bi) and 7(x 2 , b 2 ), etc. That way we obtain an expression for iV(x i, b , Y): 



• jv(xoi,bo,y) 



J ni(x i,y, bi,xi) ^7(xi,bi)^-^-d 2 &i^ + 



n 2 (xoi,y,bi,xi,b2,X2) 7(xi,bi)- — ~d b x 7(x 2 ,b 2 



d 2 xi 
2izx\' 



d 2 x 2 
2irx 2 



<Tb> 



E 



> f ( d 2 xi \ ( d 2 Xi \ 

^ J ni(xoi,y,bi,Xi, . . . ,bi,Xi) ^7(xi,bi)^— 2-d 2 6iJ . . . ^7(x i; b,)^— |-d 2 6 j; J , 



(11) 



where we put the minus sign in front of N to make it positive, since 7 is negative. Eq. (11) clarifies the 
physical meaning of iV as a total cross-section of a qq pair interacting with a nucleus. One can understand 
now the factorials in the definitions of the dipole number densities (g), (||) and ([n]): once the convolutions 
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with the propagators 7 are done then the dipoles become "identical" and we have to include the symmetry 
factors. 

In order to write down an equation for iV(xoi, bo, Y) we have to find the equations for n^s first. Following 
the techniques introduced in jl]-|} we have to differentiate the equation for the generating functional (^) with 
respect to u(x, b) putting u = 1 at the end, keeping in mind that Z(h , x i, Y, u = I) = 1. Differentiating 
formula (Q) once we obtain an equation for ni(xoi, Y, b%, xi): 



ni(x i,y,bi,xi) = <5 2 (x i -xi)27rxi 5 2 (b 1 )exp 



aCp 



dy exp 



rf 2 £ 2 2 01 2 2ni(x 2,y,bi,xi), 

;c 02 a '12 



(12) 



where, following 0], we have defined bi = b; — bo — 5X12. 

Differentiating Eq. (0) twice we obtain an equation for 712(2:01, Y, bi, xi, b 2 , x 2 ) 



n 2 (a;oi,y,bi,xi,b2,X2) 



(iy exp 



7T VP 



d 2 x 2 



'-01 



X (-19 



02-^12 



x [2 n 2 (0:02 , y , bi , xi , b 2 , x 2 ) + n% {x 2 , y , bi , xi ) m (£12 , y , b 2 , x 2 )] , 



(13) 



where b^ = b^ — bo + 5X20- Now higher order differentiation of Eq. (]?]) becomes apparent, and could be 
easily done yielding the following equation for the number density of i dipoles: 



rii(xQi,Y, bi,xi, . . . ,bi,x. 



clCf 



dy exp 



d 2 i 2 



L 02 x 12 



x [2 71,(0:02, y,bi,xi, . . . ,bj,Xj) + ^ nj(io2,y,bi,xi, . . . ,bj,Xj)nk(xx2,y,bj+x,Xj+i, ■ ■ .,bi,x;)], (14) 

where we anticipate the integration over the dipole sizes and treat the dipoles as identical objects. In 
principle Eq. ( |i!| ) should contain the permutations of the arguments of the gluon densities in the product 
on the right hand side, but for the abovementioned reason we do not write this terms explicitly. 
Multiplying formula (|l4j) by 



d 2 Xi \ / , , , d 2 x. 
W ^•••( v 7(x l ,b l )^- r ,„,. 



1 ^2 



d% 



integrating over the dipole sizes and impact parameters, and summing all such equations, i.e. summing over 
i from I to 00 in (ph one obtains the equation for 7V(xoi, bo, F) 



N(x. x,b ,Y) = -7(x i,b ) exp 



.l^lnf^i]) 

7T V P 



aC F f Y 
— — I dy exp 



7T^ 







x d x 2 



2 2 
X 02 X X2 



[2 iV(x 02 , b + -X12, y) - 7V(x 02 , b + -xi 2 , y) N(x 12 , b - -x 20 , y)}. 



(15) 



Equation (p"5|), together with equations (3) and (0), provide us with the leading logarithmic evolution of the 
F2 structure function of a nucleus including all multiple pomeron exchanges in the large- N c limit. 

Throughout the preceding calculations we never made an assumption that Q 2 is large. Of course it should 
be large enough for the perturbation theory to be applicable. The only assumption about the incoming 
photon's momentum that we made was that its light-cone component q + is large, therefore we could neglect 
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the inverse powers of q + . This is cikonal approximation, which is natural for leading ln(l/a;) calculation. 
However if the inverse power of q+ comes with an inverse power of g_, forming something like l/2q + <7_ ~ 
1/Q 2 we do not neglect these terms, therefore resumming all the inverse powers of Q 2 ("higher twist terms"). 
That way we proceed to conclude that equation ( [j"5| ) sums up in the leading logarithmic approximation 
all diagrams that include the effects of multiple pomeron exchanges, with pomeron ladders together with 
pomeron splitting vertices being incorporated in the dipole wave function. In terms of conventional (not 
"wave functional") language Eq. ( [la) resums the so called "fan" diagrams (see Fig. ||) which were summed 
up by conventional GLR equation |5|,|6j . The difference between our equation and GLR is that Eq. ( |l5| ) does 
not assume leading transverse logarithmic (large Q 2 ) approximation. 




A A A A 

FIG. 3. Multiple pomeron exchanges and splittings resummed by Eq. (^|). Each pomeron ladder interacts with a 
nucleus, which is symbolically denoted by A. 



III. DOUBLE LOGARITHMIC LIMIT 



In order to reconcile our approach with traditional results in this section we will take the large Q 2 limit 
of Eq. ([15]) and show that in this double logarithmic approximation Eq. (JlJ) reduces to the GLR equation 
||,[|. We consider a scattering of a virtual photon, characterized by large momentum scale Q, on a nucleus 
at rest characterized by the scale Aqcd- The Q 2 3> Aq CD limit implies that the dipoles produced at each 
step of the evolution in the dipole wave function must be of much greater transverse dimensions than the 
dipoles off which they were produced. Basically, since in the double logarithmic approximation the transverse 
momentum of the gluons in the dipole wave function should evolve from the large scale Q to the small scale 
Aqcd, than the transverse sizes of the dipoles should evolve from the small scale 1/Q to the large scale 
1/Aqcd- 

In the limit when the produced dipoles are much larger than the dipole by which they were produced 
(large Q 2 limit), the kernel of Eq. fllE| ) becomes 

d 2 i 2 ^^x 2 01 n ? % f) (16) 

X 02 X 12 Jx^ \ X 02J 

where X02 ~ £12 ^> x m, and the upper cutoff of the X02 integration is given by the inverse momentum scale 
characterizing the nucleus, 1/Aq CZ j. Since this integration is done in the region of large transverse sizes the 
ultraviolet cutoff p is no longer needed. One can easily see that including virtual corrections would bring in 
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the exponential factor e"^r LY in the Eqs. (Q) and (15) instead of exp - 4a ^' F In Y 



In the double 



logarithmic approximation aY \n(Q 2 / Aq CD ) > 1 and \n(Q 2 / Aq CD ) 3> 1, therefore aY < 1. That way the 
factor of e ^ Y can be neglected. The resulting limit of Eq. (|l5j) is 



rvf7^ /" y f 1 / A Qcr> ,V r 2 

N(xoi, bo, y) = -7(xoi, b ) + — / ^ / 7^2 I 2 ^( x 02, b , Y) 

n Jo Jx* (^02) 



-2V(xo2,bo,y)JV(xo2,bo,K)], 
which after differentiation with respect to y yields 



^ (X °" b0 ' r) = ^ ^ / VA " D ^ [2 N{xq2 , bo , Y) - Af(x 02 , b , y) A,(xo2, b , y)], (17) 



(x, 

where, for simplicity, we suppressed the difference in the impact parameter dependence of N on the left and 
right hand sides of Eq. (|l7). This is done in the spirit of the large cylindrical nucleus approximation. Also 
one should keep in mind that for this double logarithmic limit in the definition of N(x.Qi,h Q ,Y) given by 
Eq. ( pH ) the integration over the dipole's transverse sizes should be also done from Xq X to 1/Aq CD . 

Now we have to make a connection between JV(xoi, bo, Y) and the gluon distribution function xGa(x, Q 2 ) 
of a nucleus. iV(xoi, bo, Y) is a forward scattering amplitude of a qq pair on a nucleus and is a well-defined 
physical quantity. However there is some freedom in the definition of the gluon distribution. If one makes 
us of the general definition of the gluon distribution as a matrix element of leading twist operator, then 
an attempt to take into account higher twist operators would lead only to renormalization of their matrix 
elements (see |32| and references therein). The evolution equation for xG would be linear, with all the 
non-linear saturation effects included in the initial conditions. The goal of the GLR type of approach is to 
put these non-linear effects in the evolution equation. Therefore in the double logarithmic approach one 
usually defines the gluon distribution function through a cutoff operator product expansion, i.e., as a matrix 
element of the A^A^ operator, with Q 2 an ultraviolet cutoff imposed on the operator (see the discussion on 
pp. 442-443 of ||). In the spirit of this approach we define the gluon distribution by 

2 

AT(xoi,b , Y) = —^- x 2 01 xG A (x, l/x 2 01 ), (18) 

with the coefficient fixed by the two gluon exchange between the quark-antiquark pair and the nucleus (in 
the large N c limit). Substituting Eq. ([TJ|) into Eq. ( |l7| ) one obtains 



Differentiating the resulting equation with respect to ^(I/x^Aq^) and using xqi ~ 2/Q, which is valid in 
the double logarithmic limit, we end up with 

&xGA*,C?) _ aN , xGAix , Q2) _ <^±. [xGa {x,Q*)]\ (19) 



d\n{l/x)d\n{Q 2 /k 2 QCD ) TT >-* > s± Q 2 

which exactly corresponds to the GLR equation |5|||, with the factors matching those corresponding to 
cylindrical nucleus case in references ■ 

One has to note that the problems with the definition of the gluon distribution function outlined above bear 
no consequence on Eq. ([15]). This equation describes the evolution of AT(xoi,bo,y) in the leading ln(l/a;) 
and does not assume collinear factorization or impose transverse momentum cutoffs, therefore posing no 
problems like the mixing of operators of different twists |32| . 

One would also like to rederive the equation for xGa derived earlier in |ll],|l8|,^l|] , given by Eq (0) here. 
In the double logarithmic limit we did not reproduce those results, given by formula (48) in fll8[ ], formula 
(3.19) in pi], and, most explicitly by formula (41) in pjj. Nevertheless this equation could be obtained from 
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Eq. ( |15[) in the following peculiar limit. Consider the scattering of a virtual photon with moderately large 
virtuality Q on a nucleus composed of very heavy quarks, so that the typical momentum scale characterizing 
the nucleus given by the mass of the heavy quarks M is much larger than Q. In that case the double 
logarithmic limit would correspond to production of small dipoles in the dipole wave function. Now, using 
certain freedom we have in the definition of the gluon distribution function, we can treat Eq. (^|b) as a 
definition of xGa, which should be substituted in (^b) instead of AxG. Plugging the generalized Eq. (^jb) 
into Eq. (||) taken in the small-dipole limit we obtain Eq.(^). That way we will rederive the results of 
Jj^ , fL8|]2l| , however not quite for the same process as the one for which they were derived originally. 

Finally we note that we failed to find a limit in which Eq. (|l^) reproduces formula (39) in [^5), that was 
given by Eq. (^|) above, which is another candidate for the double logarithmic evolution of xG including all 
multiple pomeron exchanges. 



IV. CONCLUSIONS 



One of the diagrams that were not included in our analysis above is shown in Fig. ^ (a) . There a hard 
pomeron ladder splits into two and then the two ladders again fuse into one which connects to the nucleus. 
In general this diagram is of the same order in coupling and rapidity as the usual two-pomeron exchange 
diagram of Fig. ^ (b), and, therefore should be also considered. The effects of that type can play an 
important role in unitarization and saturation of the onium-onium scattering amplitude |p3[ . In the dipole 
wave function language the graph in Fig. ^ (a) can be interpreted in two ways. In one case it could be a part 
of the dipole wave function, which should be redefined to include this graph in the form of color quadrupoles. 
On the other hand the graph in Fig. [| (a) could correspond to the usual "two dipole" part of wave function 
(n.2) with both dipoles interacting with the same nucleon in the nucleus. However, for the given problem of 
scattering on a large nucleus, the graph in Fig. |J (a) is suppressed by powers of atomic number A compared 
to the graph in Fig. [| (b) for both cases. 




A A A 



(a) (b) 
FIG. 4. (a) Diagrams which are not included in our analysis, (b) A diagram which is included in Eq. 

The inclusion of the pomeron fusion effects in the dipole wave function is a difficult task @. Eq. (|) does 
not take them into account. Construction of dipole wave function which includes the diagrams shown in Fig. 
^ (a) is an interesting and important problem, which is still to be solved. One should note that inclusion of 
the graphs of the type shown in Fig. ^ (a) in the dipole wave function would result in appearance of color 
quadrupoles [j53). In the case of onium-onium center of mass scattering this wave function can be either 
considered as suppressed by factors of Nq compared to the one-dipole wave function, or, more correctly, 
one may note that its contribution to the scattering amplitude will suppressed by the factor of e^ ap ~ 1 ^ Y l 2 
||, with ap the intercept of the BFKL pomeron. However, this argument does not apply for the case when 
one of the onia is at rest Ip3l, or to our case when the nucleus is at rest. 
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The graph in Fig. |4| (a) in our case of nucleus at rest brings in suppression by powers of A. The reason for 
that is very simple: in the first diagram (Fig. |^ (a)) there is only one pomeron ladder (dipole) interacting 
with the nucleus below, whereas in the second diagram (Fig. || (b)) there are two ladders. Since each dipole 
is convoluted with its propagator through the nucleus, each ladder brings in a factor of 7(x, h)d 2 b. For large 
transverse size dipoles this factor is proportional to A 2 / 3 and for small dipoles it scales as A 1 . In any case 
one can see that the graph in Fig. |^ (a), having the same parametrical dependence on a and Y as the graph 
in Fig. [| (b), is suppressed by some power of A compared to this second graph. Therefore, one should note 
that by neglecting the effects of the graph in Fig. [| (a) we assume that we are doing the leading calculation 
in the powers of the atomic number of the nucleus A. That allows us to avoid complications associated with 
the incorporation of the diagram in Fig. ^ (a) in the evolution equation (|T5|). 

Eq. ( |T5| ) can also be derived directly from Eq. (fjj) by putting u(x i, b ) = 7(x i, b ) + 1 in it and defining 
N = 1 — Z. That way we have a method of resumming all multiple pomeron exchanges for any higher order 
corrections to the dipole kernel. If one calculates the dipole kernel, say, at the next-to-lowest order, then 
we can write down an equation for generating functional Z, similar to Eq. ([?]). Though the next-to- lowest 
order equation will in addition have cubic terms in Z on the right hand side. Then, putting u — 7 + 1 and 
Z = 1 — N one would easily obtain an equation resumming multiple pomeron exchanges in the sublcading 
logarithmic approximation. Therefore dipole model provides us with a relatively straightforward way of 
taking into account the multiple pomeron exchanges once the one-pomeron exchange contribution has been 
calculated. In other words if the dipole kernel is known at any order in the coupling constant one can easily 
generalize the resulting equation for generating functional to include the multiple pomeron exchanges on a 
nucleus. 

Finally we note that it would be interesting to try fitting the recent HERA data (li| using the evolution 
of the F2 structure function given by Eq. (|lj) . 
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